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For any mean value of a cartesian component of a spin 
vector we identify the smallest possible uncertainty in any of 
the orthogonal components. The corresponding states are op- 
timal for spectroscopy and atomic clocks. We show that the 
results for different spin J can be used to identify entangle- 
ment and to quantify the depth of entanglement in systems 
with many particles. With the procedure developed in this 
letter, collective spin measurements on an ensemble of par- 
ticles can be used as an experimental proof of multi-particle 
entanglement. 

The commutator relation for angular momentum op- 
erators leads to Heisenberg's uncertainty relation for the 
Cartesian components 



AJ, •AJy> |(J,)|/2. 



(1) 



Without violating Heisenberg's uncertainty relation, it 
is possible to redistribute the uncertainty unevenly be- 
tween Jx and Jy, so that a measurement of either 
or Jy become more precise than the standard quantum 
limit y/\{Jz) 1/2. States with this property are called spin 
squeezed states in analogy with the squeezed states of a 
harmonic oscillator. 

A two-level atom can be represented as a spin 1/2 sys- 
tem, and in experiments on a large number of atoms N 
which all start out in the same initial state and which 
are all subject to the same Hamiltonian, one can conve- 
niently express the collective observables of the gas by 
means of an effective spin J = A^/2, so that, e.g., the 
difference in number of atoms populating the two in- 
ternal states is given by the J^-component. The state 
with all atoms in the "spin up" internal state, is equiva- 
lent to the \Jz = J) eigenstate of the macroscopic spin. 
If one measures the x component of the spin of a sin- 
gle atom, it is projected onto the internal superposition 
states (I I) ± I 1))/^/2 with equal probability. The value 
of the total Jx is given by the difference in the number of 
atoms in the two states, and it fluctuates binomially with 
a variance J/2, which matches precisely the equality sign 
in (|lj) with AJx = AJy. 

The states and the amount of spin squeezing produced 
by applying Hamiltonians and ~ Jy to an initial 
\Jz = J) state have been studied 0|, and the squeezed 
states which satisfy the equality sign in (0), the so-called 
minimum- uncertainty- product states, have been identi- 
fied PI |4| . Interaction of atoms with broadband squeezed 
light MM, is an experimentally verified means to produce 
spin squeezing 0. Spin squeezing based on ideas from 
quantum computing was recently suggested Q, and re- 
cent ideas for neutral atom spin squeezing based on QND 



detection of the atomic spin state JP| and on the colli- 
sional interactions between atoms @jll) give reason to 
believe that sizable spin squeezing may be much easier 
to achieve than optical squeezing. 

In Ramsey type spectroscopy on a collection of two- 
level atoms, a signal proportional to the length of the 
mean collective spin pointing, say, along the z-axis is 
recorded and the noise is given by the uncertainty of one 
of the orthogonal components. Wineland et al. have 
shown that the frequency resolution in spectroscopy 
on N two-level atoms contains the factor 



2JAJx 



(2) 



which is reduced by spin squeezing. In this way, spin 
squeezing becomes an important ingredient in high pre- 
cision spectroscopy and in atomic clocks, which are at 
present limited precisely by the fundamental spin noise 
JlSf . Furthermore, spin squeezing is an important ingre- 
dient in quantum information, because the ensuing quan- 
tum entanglement leads to possibilities, e.g. for atomic 
teleportation [l^ . 

In the derivation of Eq. it is assumed that no other 
sources of noise are present. The states which minimize 
the quantity ^ are obtained in the limit (J^), AJ^. 0, 
where any other source of noise will, however, deteriorate 
the spectroscopic resolution. In practice, the ideal states 
for spectroscopy are therefore states which minimize the 
noise AJx for a given not too small value of (J^). In 
the present Letter we identify this minimum, i.e., we op- 
timize the signal to noise ratio by identifying the states 
with minimum quantum noise for a given value of the sig- 
nal. For photons a similar analysis has been performed in 
Ref. [ p^ . Having identified the minimum of AJx we use 
this information to derive an experimental criterion for 
entanglement. A measurement of two collective operators 
for a collection of atoms produces a sufficient criterion 
for entanglement which can even quantify the "depth" 
of entanglement, i.e., the minimum number of particles 
forming multi-particle entangled states in the sample. 

To get a lower limit on AJx as a function of (Jz), one 
can use the inequality {Jx) + {Jy) + {Jz)'^ < JiJ+}) which 
together with Heisenberg's uncertainty relation (n|) yields 
the limit 



{AJx)''>^ J{J+l)-{Jz? 



^^{J{J+l)-{Jz)^f-{JzY 



(3) 



This does not present a tight minimum for Var(Ja:), but 
for large J and ( Jz) ~ J it is close to the actual minimum 
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found by the numerical approaches discussed below. For 
low values of ( Jz) it differs by a factor of two. The pre- 
cise analysis of the minimum becomes quite different for 
integer spins and for half-integer spins, and we shall deal 
with them separately: 

For integer spins our calculations show that the state 
which minimizes Var(Ja:) for a given {Jz) has vanishing 
{Jx) and (Jy), so that it is also a minimum of (J|). 
Accordingly, these states can be found by minimizing 
p^{Jz) + {Jx)^ where /x is a Lagrange multiplier, ensur- 
ing the value of {Jz)- For J- values up to several hun- 
dred, it is straightforward to numerically determine the 
minimum, by determining the smallest eigenvalue of the 
operator ^ -I- for a wide range of values of ^. By de- 
termining {Jz) and {Jx) in the corresponding cigcnstatcs 
one finds exactly the minimum value of Var( Ja;)=( J^) 
for the particular value of {Jz)- The results for differ- 
ent values of J are shown in Fig. ||. For J = 1 it is 
possible to diagonalize /i -I- analytically, and we get 
Var( J,)„„„ = (1 - ^1 - ( Jz)2)/2. 
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FIG. 1. Maximal squeezing for different values of J. The 
curves starting at the origin represent the minimum obtain- 
able variance as a function of the mean spin. Starting from 
above, the curves represent J = 1/2, 1, 3/2, 2, 3, 4, 5, 10. The 
dotted curve for J = 1/2 is the limit identified in Ref. 0. 
The full curves are obtained by diagonalization of the opera- 
tor /xJz -I- . The dashed curve represents the position of a 
bifurcation in the solution for half-integer spins. To the right 
of this curve the diagonalization may be applied. To the left 
of the curve the minimum is found by a variational calculation 
(dash-dotted curve for J = 3/2). 

For half-integer spins, it is not true that the state min- 
imizing (J^) also minimizes Var( J^,) at a given value of 
(Jz). The reason is that for half- integer spins, the op- 
erator has eigenvalues ^, |,..., and its expectation 
value will thus always exceed \. The variance of J^, 
however, can come arbitrarily close to zero, if the sys- 
tem approaches a Jx eigenstate. Consider for instance a 
J = \ particle, where all (pure) states can be obtained 
as a simple rotation of the spin up state. In this case the 
components perpendicular to the mean spin are never 



squeezed; their variances are both i. But if we compute 
the variance of Jx and the mean value of Jz, one finds 
that they obey the relation, Var(J2;)mm = {Jz)"^ -, where 
both sides approach zero when the state approaches ei- 
ther of the two Jx = ^\ eigenstates. In that case, of 
course, the mean spin also has a component along the x 
direction. The state is spin squeezed in the sense of the 
relation (|^), but not in the sense where one deals explic- 
itly with a spin component perpendicular to the mean 
spin vector. 

For large half-integer J it is more difficult to find the 
most squeezed states. The reason is that the problem 
cannot be formulated as a linear quantum mechanics 
problem like the diagonalization of an operator contain- 
ing a Lagrange multiplier term, which we used for integer 
spins. To compute a variance, we have to determine the 
square of a mean value which is an expression to fourth 
order in wave function amplitudes. It is easy, however, to 
implement a Monte Carlo variational calculation which 
minimizes /i(Jz) + Var(Ja;), by randomly modifying the 
amplitudes of a state vector as long as the variational 
functional is reduced. Like above, the Lagrange multi- 
plier term is used to adjust the mean value of Jz, so that 
for each value of [l the identified state vector minimizes 
Var(J2;) for the given value of (Jz). When applied to 
larger half-integer values of J this method shows that 
large values of (Jz) are accompanied by vanishing mean 
values of Jx and Jy , and the solution thus coincides with 
the one found by the diagonalization method. But, for 
a critical value of (Jz), the solution bifurcates, and two 
states with opposite nonvanishing mean values of Jx have 
the smallest variance. See Fig. |[ These states approach 
the two Jx = ±5 in the limit of small (Jz). Due to 
the noise in the simulation, the Monte-Carlo method is 
not efficient for a precise determination of the critical 
point of the bifurcation for large values of J. Before the 
bifurcation the state is the eigenstate corresponding to 
the lowest eigenvalue of the operator /iJz -f J^, and after 
the bifurcation the state is a superposition of the differ- 
ent eigenvectors with amplitudes on states with higher 
eigenvalues. We can therefore determine the position of 
the bifurcation from the properties of the eigenvectors, 
and for different values of J < 100, we find that the bi- 
furcation happens in the interval 

0.83 < < 0.88, (4) 

except for the special case J = 1/2 where (Jz) = J at 
the bifurcation. If we do not break the iJ^; symmetry 
the variance from this point flattens out to the value -j, 
but in either of the states with the broken symmetry, 
the variance decreases towards zero when smaller values 
of (Jz) and non-vanishing values of {J^ are considered. 
The position of the bifurcation is plotted in Fig. |l|. To 
the right of the dashed line the minimum may be found by 
diagonalization. To the left of the curve the variational 
approach has to be applied for half-integer spins. 
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FIG. 2. Bifurcation in the solution for J = 3/2. The 
points represent the mean value of the spin in the maxi- 
mally squeezed states. The points are obtained by a Monte 
Carlo variational calculation which minimizes /x( Jz)-|-Var( Jj;). 
Above (Jz) ~ 1.32 there is a unique solution with {Jx) = 0. 
Below (Jz) ~ 1.32 the solution bifurcates, and {Jx) ap- 
proaches ±1/2. 

It is the non-linearity of the problem that leads to the 
bifurcation and symmetry breaking of the solution. Clas- 
sical approximations to many body quantum problems 
often show such bifurcations associated with phase tran- 
sitions in the problem, e.g., lasing. It is interesting that a 
similar phenomenon appears here, in the study of a sin- 
gle quantum system in a (very) low-dimensional Hilbert 
space. We emphazise that we are not discussing an exten- 
sion of quantum theory to include non-linear terms, we 
are simply identifying the quantum states that minimize 
a variance, and this is a non-linear problem. 

Since we have identified the maximally squeezed states 
as eigenstates of the operator ^Jz + Jj, even without 
having explicit expressions for these eigenstates, we can 
devise a method to produce them. This method only 
works for integer spin, and for half-integer spin which 
are squeezed to values of the mean spin exceeding the 
value at the bifurcation (Q). The system is first pre- 
pared in the | Jj = J) eigenstate, and one switches on a 
Hamiltonian H{t) = ujJz +x{t)Jx' where increases 
very slowly from the value zero and where w < 0. If the 
state of the spin follows this Hamiltonian adiabatically, 
it evolves through the minimum energy eigenstates of the 
instantaneous H{t), which is precisely the family of states 
identified by the above diagonalization procedure. The 
adiabatic process may be difficult to perform in physical 
systems of interest, and for practical purposes it is rele- 
vant to point out that the straightforward application of a 
Hamiltonian H = J^, also leads to spin squeezing [Q and 
in the regime with large |(Jz)|, the squeezing resulting 
from this Hamiltonian is actually close to the optimum. 
The Hamiltonian H = J^ — Jy, also discussed in Ref . [Q , 
leads to similar squeezing for large KJz)], and it follows 
the optimum for a larger range of parameters. 

We have identified the minimum variance of given 
the value of (Jz). Any measurement of these two quan- 
tities can be plotted as a point in Fig. |l|, and this point 
must lie on or above the curve for the relevant J. We 
note that the curves depend on J, and in the chosen 
units, large spins can be more squeezed than small spins. 



This implies that the collective spin variables J = '}2iiJi 
for several spin J particles can be more squeezed than the 
individual spins themselves. We will now show that this 
requires the state of the spins to be an entangled state. 
It is already known jl^, that for spin 1/2 particles, re- 
duction of the parameter ^ below unity for the collective 
spin implies entanglement of the spins. We generalize 
this property to arbitrary spins. 

A separable (non-entangled) state of N spin J parti- 
cles is defined as a weighted sum of products of density 
matrices with positive weights pk p6|,p^ : 



E(fc) ^ (k) ^ (k) 



(5) 



where p'^^ is the density matrix of the i*^ particle in the 
k*-^ term of the weighted sum. The variance of Jx in such 
a state obeys the inequality 

N 



k i=l 
N 



(6) 



k i=l 



where the function Fj{-) is the minimum variance of Jx 
divided by J for the spin J particle, i.e. the curves plot- 
ted in Fig. I, and (Jz)f ^ is the mean value of Jz of the 
i*'^ particle in the fc*'* realization in the sum (H). 

As it appears from Fig. ^ all the curves Fj{-) are con- 
vex. We can prove this property for integer spins, and for 
half integer spins in the range of large \ {Jz)\, by consid- 
ering the production of the states by adiabatic passage 
from the \Jz = J) eigenstate. The positive factor in front 
of the component in the Hamiltonian ujJz +x{'t)Jx 
gradually increased, and the rate of change of ( Jj) and of 
( Jz) at any given time are given by Ehrenfest's theorem: 



(7) 



The mean values on the right hand side should be evalu- 
ated in the maximally spin squeezed state, i.e., they are 
not known explicitly. But, we observe that they contain 
the same operator, and the ratio between the two rates of 
changes is therefore simply —oj/x{t). This implies that 
along the family of maximally squeezed states, the rela- 
tive change of (Jj) and (Jz), i.e., the slope of the curve 
Fj(-), is monotonically increasing (since x(t) is an in- 
creasing function of time and u < 0). It follows that the 
second derivative of the function Fj{-) is positive, i.e., 
the function is convex. 

From the convexity follows that the functions Fj{-) 
obey Jensen's inequality, which states that any linear 
combination of Fj{ai)'s with positive coefficients is larger 
than or equal to the function Fj evaluated on the linear 
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combination of the arguments. It therefore follows that 
in any separable state 

Var(j7.) > Y.PkNJFj f](J.)f M 

>ivjf.(i:»i7i:(A>!") (8) 

This relation is a main result of this paper. In an ex- 
periment with a collection of N spin J particles, it is 
possible to measure the collective Jz and J^^ and to de- 
termine their mean value and variance. If the data-point 
{{Jz)lNJ, Ya,v{Jx)lNJ) lies below the relevant curve in 
Fig. ^, the systems are not in a separable state, i.e., they 
are experimentally proven to be in an entangled state. 

The extent to which the measured data point falls be- 
low the curve in the plot is a measure of the degree of 
entanglement. A quantitative measure of entanglement 
in a multi-particle system is the number of elements that 
must at least have gone together in entangled states. We 
define a k particle entangled state to be a state of N 
particles which cannot be decomposed into a convex sum 
of products of density matrices with all density matrices 
involving less than k particles: at least one of the terms 
is a /c particle entangled density matrix. If, for example, 
N spin i particles form N/2 entangled pairs, the degree 
of macroscopic spin squeezing of the system is limited by 
the inequality (||) with J = 1 and N replaced by N/2. 
If the measured macroscopic {Jz) and Yai{Jx) also lie 
below the corresponding J = 1 curve, the measurement 
unambiguously implies that the systems are entangled in 
larger than binary ensembles. The size of these ensem- 
bles is a measure of the "depth" of entanglement, which 
can be determined experimentally This criterion may 
be compared to the one used in [|8| where the fidelity of 
production of a maximally entangled A'^-particle states is 
used as a proof of A^-particle entanglement. 

As a final point we demonstrate how our procedure can 
be applied to identify substantial multi-particle entangle- 
ment in recent theoretical proposals for spin squeezing. 
In Ref. PI it is predicted that appreciable spin squeez- 
ing of atoms can be obtained in a two-component Bose- 
Einstein condensate. In the calculation, a reduction of 
N&t{Jx) by a factor of 1000 is found for a reduction of 
{Jz) of only 1% with 10^ atoms in the condensate. Using 
Eq. these numbers imply a depth of entanglement of 
~ 2 • 10**. In ion traps it has been shown that it is possible 
to implement a Hamiltonian by applying bichromatic 
light to all ions in the trap This Hamiltonian can 
be used to create a maximally entangled state of all the 
ions. If the decoherence in the trap is such that one 
cannot produce a maximally entangled state, a different 
strategy could be to apply the light for a short time so 
that squeezing is produced. For small times the squeezing 



obtained from the Hamiltonian is close to the optimal 
curves in Fig. [|, and our theory identifies a depth of en- 
tanglement close to the total number of ions in the trap. 
In this way one could produce and verify the production 
of an entangled state of many ions. 

We have considered squeezing and entanglement re- 
lated to collective vector operators Jz and Jx- We em- 
phazise that the collective spin components of multi- 
particle atomic system are readily available by standard 
spectroscopic methods, which require no access to the 
individual components. Given the large interest in spin 
squeezing, a criterion of entanglement based on this prop- 
erty is an important tool. Recall, however, that systems 
may well be entangled without being spin squeezed: The 
spin squeezing measurement provides a sufficient crite- 
rion for the depth of entanglement, not a necessary one. 
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